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ABSTRACT 

We examine the effects of environment on the dynamical structure of sateUite systems 
based on the Millennium-II Simulation. Satellite halos are defined as sub-halos within 
the virial radius of a host halo. The satellite sample is restricted to those sub-halos 
which showed a maximum circular velocity above 30 kms~^ at the time of accretion. 
Host halo masses range from 10^^ to 10^* /i~^Mq. We compute the satellites' average 
accretion redshift, Zacc, velocity dispersion, a, and velocity anisotropy parameter, /3, 
utilising stacked satellite samples of equal mass hosts at similar background densities. 
The main results are: (1) On average satellites within hosts in high density environ- 
ments are accreted earlier {Az w 0.1) compared to their counterparts at low densities. 
For host masses above 5 x 10^'^ h~^MQ this trend weakens and may reverse for higher 
host masses; (2) The velocity dispersion of satellites in low density environments fol- 
lows that of the host, i.e. no velocity bias is observed for host halos at low densities 
independent of host mass. However, for low mass hosts in high density environments 
the velocity dispersion of the satellites can be up to ^ 30% larger than that of the 
host halo, i.e. the satellites are dynamically hotter than their host halos. (3) The 
anisotropy parameter depends on host mass and environment. Satellites of massive 
hosts show more radially biased velocity distributions. Moreover in low density envi- 
ronments satellites have more radially biased velocities (A/3 > 0.1) compared to their 
counterparts in high density environments. We believe that our approach allows to 
predict a similar behaviour for observed satellite galaxy systems. 

Key words: methods: N-body simulations — methods: numerical — dark matter — 
galaxies: haloes — galaxies: clusters: general 



1 INTRODUCTION 

The dependence of halo statistics on a second parame- 
ter in addition to mass is now generally referred to as 
assembly bias. Simple ext ensions to the Press-Schechter 
and e xcursion set models ( Press & Schcchtcr 1974 |; iKaiseJ 
Bond et al.lll99ll: [Cole fc Kaiser .1989; Lacev fc Cold 
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Mo fc White! 19961 ) predict the clu stering of halos to 
depend on their mass alone. However, iGao et al] (|2005l ) 
and various subsequent studies showed that clustering also 
depends on other halo properties, for example, formation 
time, concentration, substructure content, spin and shape 
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demonstrated that excursion set models predict that halos 
in denser environments do form later, independent of halo 
mass. At the high mass e nd this agrees with findings from 
N-bo dy simulations (e.g., IWechsler et al.l I2OO2I : Ijing et al] 
|2003), but is opposite to the behaviour observed at the low 
mass end. Based on st atistics of the peak s within Gaussian 
random fluctuations, iDalal et al.l (|2008l ) argued that the 
behaviour for low mass halos can be understood if cessation 
of mass accretion is t aken into account. A similar argument 
has been proposed bv lHahn et afl (|2009l ). 



Several other studies have investigated the dependence 
of halo formation times or, similarly, merger ra tes on envi- 
ronm e nt (e.g..,Gottlober et al.,i2001i.i2002: Shcth fc TormenI 
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I2OO4I : iFakhouri fc Mai |2009| . I2OI0I : iHahn et all |2007| ). Al 

though slightly different density estimators are employed, 
such as over density in a sphere or mark correlation func- 
tions, it is generally agreed upon that halos less massive 
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than ~ 10^'' H^^Mq, which reside in high density regions, 
form earlier compared to those with equal mass but located 
in less dense regions. 

In a preceding study we discussed the dependence of 
the dy namical structure of dark m atter halos on environ- 
ment (jFahenbacher fc WhiQ I2OI0I '). We have found that 
the velocity dispersion of dark matter halos residing in less 
dense environments shows a more radially biased velocity 
structure compared to halos of the same mass in denser 
environments. We now extend this study to the dynamics 
of satellites within host halos. Satellites are selected on 
the basis that our findings can be applied to observed 
satellite galaxy systems. The effects of environment on 
the dynamical structure of satellites may have important 
implications for the accretion processes of the host galaxies. 

The outline of the paper is as follows. In §[2] we intro- 
duce the simulation and the halo finding procedure. We then 
explain how the environment of host halos is determined and 
describe the stacking procedure to compute average proper- 
ties of satellite populations. § [3] presents the results: host 
halo properties as a function of environment are shown in 
§ 13.11 main results are portrayed in § 13.21 where we discuss 
the accretion history and dynamical structure of satellite 
systems. In ij |3.3l we investigate resolution effects. A conclu- 
sion is given in § 2] 



2 METHODOLOGY 

In this section we describe the tools for our analysis. A short 
account on the Millennium-II Simulation is followed by the 
discussion of the halo finding procedure. We also explain 
the determination of the background density and the stack- 
ing procedure used to derive averaged properties of satellite 
systems in different environments. 

2.1 Simulation 

The Millenniu m-II Simulation (|Bovlan-Kolchin et ai]|2009l : 
ISpringe]|2005l ) adopted concordance values for the param- 
eters of a flat A cold dark matter (ACDM) cosmological 
model, Qn-i = 0.25 and ilb ~ 0.045 for the current matter 
and baryon densities, h — 0.73 for the present dimensionless 
value of the Hubble constant, ag = 0.9 for the rms linear 
mass fluctuation in a sphere of radius 8 /i~^Mpc extrapo- 
lated to 2: = 0, and n = 1 for the slope of the primordial 
fluctuation spectrum. The simulation followed 2160'^ dark 
matter particles from z — 127 to the present day within a 
cubic region 100 /i~^Mpc on a side resulting in individual 
particle masses of 6.9 x 10^ h~^MQ. The gravitational force 
had a Plummer-equi valent comoving softening of 1 /i~^kpc. 
We refer readers to ^ Bovlan-Kolchin et al.l (|2009h for more 
detailed description of the simulation. 

2.2 Halo sample and accretion times 

The halos are found by a two-step procedure. In the 
flrst step all collapsed halos with at least 20 particles are 
identifled using a friends-of-friends (fof) group-flnder with 
linking parameter b — 0.2. These objects will be referred 



to as fof-halos. Then post-processing with the substructure 
algorithm SUBFIND (|Springel et allbOOll l subdivides each 
fof-halo into a set of self-bound sub-halos. Sub-halos with 
lest than 20 particles are discarded. Here we consider the 
most prominent sub-halo within each fof-halo as host halo. 
To characterise a host halo we use: the maximal circular 
velocity, Vmax, which peaks at a radius, -Rmax; the concen- 
tration, were iZvir denotes the virial 
radius, i.e., the radius of a sphere centred at the potential 
minimum of the host which comprises a mean density of 94 
times t he critical value . For a Navarro-Frenk- White density 
profile (|Navarro et al.l Il997l 'l. J?max is a factor of ~ 2.1 
larger than the scale radius Rs, which has conventionally 
been used to compute the concentration; and the three 
dimensional velocity dispersion, (T3dh, wich is based on the 
velocities of all host halo particles after subtracting their 
common bulk velocity from each of them. The velocity 
dispersion of each individual halo, itsdh, is scaled its virial 
velocity Vvir ~ (GMvir/fivir)^, wherc A/vir is the virial 
mass and G is the gravitational constant. 

Besides the host halo, all other sub-halos with more 
than 100 particles and located within the virial radius of 
the host are referred to as satellite halos. We only take into 
account sub-halos within the virial radius of the host halo to 
prevent the contribution of sub-halos within falsely linked 
'dumbbell' shaped fof halos. The lower limit of 100 particles 
is imposed to guarantee a reliable determination of the satel- 
lites' velocities. The accretion redshift of a satellite halo onto 
the host halo is determined as the redshift at which it dis- 
continues to be the most prominent sub-halo in its own fof- 
halo. The time elapsing between two consecutive snapshots 
is roughly 0.3 Gyr which determines the accuracy the accre- 
tion times used here. We define the redshift of the snapshot 
just prior to accretion as accretion redshift. This typically 
corresponds to the time when the halo acquires the maxi- 
mum mass during the course of its evolution (cf.. lGuo et al.l 
I2OI0I ). With Knax.acc we deuotc the maxim um circular ve- 
locity of the halo at that time. Following iKravtsov et al.l 
((2OOJ) we include only those satellite halos in our analy- 
sis with T/max,acc > 30kms~^. This restriction results in a 
mean of few 1000 simulation particles per arriving satellite. 
Halos with lower masses presumably contain only negligible 
amounts of stars. We abandon this requirement only for a 
resolution study in § 13.31 



2.3 The background density field 

The background density field is computed based on all 
sub-halos with maximum circular velocities, Vmax, be- 
tween 200 kms~^ and 300 kms~^. As Figure [T] illustrates, 
this velocity range corresponds to halo masses close to 
Mt = 6.15 X 10^^ /i~^M0 which is the typical collapse 
mass at 2 = for the given cosmology. Halos in this 
mass range provide a relatively unbiased sampling of the 
overall density field. We denote these halos as l/«-halos. 
To facilitate the comparison with density estimates based 
on observed galaxies we include all sub-halos within the 
given velocity range irrespective of whether they are host 
or satellite halos. In total there are 2013 Vt halos resulting 
in a sampling of the density field based on a point set with 
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Figure 1. Average host halo properties, maximum circular velocity (Knax), concentration (cmax) and 3D velocity dispersion ((T3oh)i as 
a function of halo mass. Red (solid) and blue (dot-dashed) lines display results based on halos in the upper and lov^er 33% tails of their 
background density distributions. Halo concentrations, Cmax, are defined as the ratio of the virial mass to the radius of the maximum 
circular velocity. Before averaging the velocity dispersion of each individual halo is scaled by its virial velocity, Vvir = (CA/vir/iJvir)' . 
The shaded regions indicate the Icr confidence intervals based on a bootstrap resampling. 



a mean distance of slightly less than 8 h ^Mpc. 

The background density for any given halo is 

determined by the seven nearest Vt halos, each of 

them smoothed with a sm oothing kernel of the form 
l|Monaghan fc Lattanziolll985l ') 
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As in lSpringel et al] (|200lh we define the smoothing kernel 
on the interval [0, h\ and not on [0, 2h\ as frequently done 
in the literature. The background density at a given halo 
location is the sum of the contributions of the seven near- 
est neighbours. Using this recipe, a background density is 
assigned to each host halo which we also address as environ- 
ment of the host. In the subsequent analysis the host halos 
are ordered according to their background density and av- 
erage properties of the upper and lower 33% tails are deter- 
mined separately. We also investigated the outcome based 
on the upper and lower 20% tails and found very similar re- 
sults. To improve statistics we choose the larger background 
density intervals. 

2.4 The stacking procedure 

Besides deriving the average accretion redshift, Zacc, of 
the satellite populations as a function of host mass and 
environment, we aim to study differences in their average 
dynamical properties. To achieve reasonable statistical 
significance individual host halos with similar masses 
and background densities are stacked. For that purpose 
we subtract the bulk velocity of the host halo from the 
individual satellite velocities and scale them by the virial 
velocity of the host, Vvir- The scaling is done to compensate 
for the host mass variation within the individual mass bins. 

The scaled velocities of satellites belonging to host ha- 
los of a given mass at a given background density are used 
to compute the mean radial velocity, Urad, the three dimen- 
sional velocity dispersion, ctsd , its radial and tangential com- 
ponents, CTrad and atan, and based on these quantities the ve- 
locity anisotropy parameter, /3 = 1 — 0.5 (f7"tan/c"rad)- In addi- 



tion we compute the average number of satellites, A'sat , as a 
function of host mass and environment as well as the average 
ratio between the actual maximal circular velocity and that 
one at the time of the satellites' accretion, Knax/'Knax,acc. 



3 RESULTS 

Before we discuss the environment dependence of accretion 
times and dynamical structure of satellite systems we re- 
view the impact of environment on the host halo properties 
themselves. Resolution effects are discussed in the final para- 
graph. 



3.1 Host halo properties as a function of 
environment 

The left panel of Figure [1] shows the average maximum cir- 
cular velocity, Knax, of host halos as a function of mass. Red 
(solid) and blue (dot-dashed) lines, as with all the remain- 
ing figures, show the result for the halos within the upper 
and lower 33% tails of the background density distributions. 
Subsequently, we will address the two samples as high and 
low density samples. The figure indicates that the average 
values of Knax are independent of environment. 

The middle panel of Figure[T]displays the concentration, 
Cmax, as a function of host mass and environment. In agree- 
ment with previous studies we find low mass halos with high 
background densities to be more concentrated compared to 
their low density counterparts of the same mass. Again, this 
behaviour is reversed for halo masses abo ve ^ l O" fe-^Mp 
(cf., IWechsler et al. I l2006l: iGao fc White! l2007l: Ijing et al 



goo?': 'Maccio et al.''2007': 'Wetzel et al.ll2007l : lAngulo et~ 

2008; Faltcnbachcr & White 2010). 

The right hand panel of Figure [T] depicts the three di- 
mensional velocity dispersion of the host halos, ctsdh, as a 
function of mass and environment. Low mass halos in high 
density regions have higher velocity dispersions compared to 
their counterparts in low density regions. This behaviour is 
reversed for higher masses which is similar to what is seen 
for the concentrations. The analog behaviour of concentra- 
tion and velocity dispersion is in a greement with predictions 
based on the Jeans equation (e.g.. lFaltenbacher fc Mathews! 



4 Faltenbacher 




Figure 2. Average sub-halo halo properties, redshift of accretion (zacc), number of sub-halos per host halo (A^sat), the ratio of the current 
maximum circular velocity and its values at the time of accretion (Vmax/Vmax.acc), the radial velocity (frad)i the three dimensional velocity 
dispersion (0-30) and the velocity anisotropy /3 as a function of the host halo mass. Red (solid) and blue (dot-dashed) lines display results 
based on host halos in the upper and lower 33% tails of their background density distributions. (By analogy, the upper and lower lines 
of the pairs of dashed and dotted lines correspond to results based on host halos in the upper and lower 33% tails of their background 
density distributions.) Before averaging the velocity of each individual sub-halo is scaled by the virial velocity of its host halo. The shaded 
regions indicate the la confidence intervals based on a bootstrap resampling. 



I2OO7I ). However, we notice the crossing takes place at some 
lower masses than seen for the concentrations. 



3.2 Accretion history and dynamical structure of 
satellite systems 

The upper left panel of Figure [5] displays the average ac- 
cretion redshift of satellites within host halos at high (red, 
solid line) and low (blue, dot-dashed line) densities as a 
function of mass. For hosts below the typical collapse mass, 
M* — 6.15 X 10^^ ft.~^Mo, satellites in high density regions 
are accreted earlier compared to satellites in low density 
regions. For host masses above M* the difference becomes 
marginal or may even reverse. The dependence of satellite 
accretion times on environment is simila r to the observe d 
assembly bias for halos themselves (e.g.. [Cao et al.ll2005h . 
A similar correlation between host halo mass and satel- 
lite accretion reds hift has also been reported in a study by 
iLagos et all (|2009l ). This behaviour can be explained by the 
fact that as the host mass increases the mean mass of the 
accreted satellites increases as well and therefore makes the 
satellite population as a whole more resistant against tidal 
forces. 

At the first glance our fin dings seem the contradict re- 
sults by Fakhouri fc Mai (|2009l ). They find a modest increase 
of the merger rates in high density environments at present 
time. However, the measurement of the mean accretion red- 
shift of surviving satellites is a convolution of the merger 
rate (over an extended period of time) and the survival rate 
which may depend on environment as the velocity struc- 



ture does. In particular we will see below that satellites in 
high density environments show a more tangentially biased 
velocity structure which prevents rapid destruction by the 
hosts' tidal field as it is the case for more radial orbits. We 
also checked, whether the satellites' masses at the time of 
accretion are dependent on environment. It turned out that 
this is not the case. Therefore, a dependence of the satel- 
lite masses on environment can be excluded as explanation 
of the accretion time differences of the surviving satellite 
population. 

The mean number of satellites for the given host halo 
masses does not depend on environment as the overlapping 
graphs in the upper middle panel of Figure [2] indicate. As 
shown in the upper right panel the average decrease of the 
maximum circular velocity is about 15% which is slightly re- 
duced for low mass hosts. Consequently, the central profiles 
of satellite halos remain fairly intact on their orbits within 
the potential well of the host. This behaviour is i n a agree- 
ment with the study by Kazantzidis et al] l|2004) , see also 
iBoylan-Kolchin et al.l (|2009l ). 

The lower left panel of Figure [2] depicts the mean radial 
velocity of satellites as a function of host mass and environ- 
ment. All velocities are scaled by the corresponding virial 
velocity. A radial velocity close to zero indicates that on av- 
erage the satellite distribution is static which is observed for 
host halos above 10^^ /i~^Mq. Negative values correspond 
to an net infall or some other excess of inward moving satel- 
lites. 

The lower middle panel in Figure [2] displays the 3D 
(solid and dot-dashed lines) velocity dispersion of satellites 
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Figure 3. Satellite velocity dispersion scaled by the host halo 
velocity dispersion, 530, as a function of host halo mass. Except 
for the scaling this figure repeats the solid and dot-dashed lines 
in the velocity dispersion panel of Fig. (2] 



and their radial (dotted lines) and the tangential (dashed 
lines) as a function of host mass and environment. All satel- 
lite velocities are computed relative to the bulk velocity of 
the host halo and scaled by the corresponding virial velocity 
before averaging. We find a remarkable rise in the disper- 
sions of the satellites in low mass hosts at high densities. 
These satellites show velocities 30% larger than their coun- 
terparts in low density environments. This behaviour is dis- 
cussed more thoroughly below in the context of Figure O 

The lower right panel of Figure [5] presents the 
anisotropy parameter of satellite galaxies as a function of 
host mass and environment. We find a profound difference in 
the dynamical structure of satellite systems in high density 
environments as compared to their counterparts at low back- 
ground densities (A/3 ^ 0.1). Independent of host mass, the 
satellite velocities in low density environments are more ra- 
dially biased than those in high density regions. Besides the 
dependence on environment there is also an obvious trend 
with host mass. Independent of environment, the velocity 
structure of satellites is more radially biased within more 
massive hosts. 

Host halos in a low density environments dominate 
the ambient gravitational field more strongly than their 
equal mass counterparts at high densities. The gravitational 
field lines point radially towards the host and so does the 
acceleration exerted on future satellites. The dynamical 
structure of the satellites within the host halo is a result of 
the fairly radial inflow. In contrast, the gravitational field 
in high density regions is more complex. Before satellites 
are accreted onto a host they also experience non radial 
accelerations from other massive haloes nearby. As a result 
the velocities of satellites show a larger non-radial compo- 
nent. The dependence of P on mass may be interpreted in 
the same way. The more massive the host the more radial 
the gravitational acceleration of future satellites. Which 
results in the more radially biased velocity dispersions of 
the satellites within high mass hosts. 

Figure [3] further explores the large velocity dispersion 
of satellites in low mass hosts at high background densities. 



Here we display the three-dimensional velocity dispersion 
of the stacked satellite populations scaled by the measured 
three dimensional velocity dispersions of the host halos and 
not, as in Figure [2] by the virial velocities. Thus, the ve- 
locity dispersion of the satellites are directly compared to 
that of the host. We find that satellites in host halos at 
low densities and satellites in high mass hosts (independent 
of environment) don't show velo city bias. This i s in a gree- 
ment with resu lts in iFaltenbacher fc Diemandl l|2006l l and 
iLau et al.l(|2009l '). The velocity dispersion of satellites in high 
density environments deviates from this behaviour. For host 
masses below ~ 10^^ /i~^M0 the satellites are hotter than 
the overall dark matter component. 

This effect can not be explained by more recent accre- 
tion times since, in this case, the low den sity satellites should 
show even larger di spersions. Recently, IWang et al.l (|2007l ) 
and Fakh ouri fc Ma (2009) reported on dynamically hooter 
ambient environments for halos in high density regions. In 
addition, low mass hosts in dense environments are located 
next to more massive halos. These halos expel a s ubstantial 
fraction o f their sub -halos out to large radii (e.g.. lGill et al.l 
|2004; Lu dlow et al. 1 12009 ). Thus, satellites which merge onto 
or pass through the low mass host systems are hotter com- 
pared to their counterparts at low densities. It seems plau- 
sible that some fraction of satellites within low mass hosts 
at high densities do originally not belong to the Lagrangean 
volume of the host halo instead they are accelerated by more 
massive structures nearby. These satellites are dynamically 
hotter then the host itself. In the following discussion we 
address these satellites as interlopers. 



3.3 Abandoning the mass barrier for entering 
satellites 

Figure |4] iterates the analysis presented in Figure [2] with the 
only difference that the restriction on the maximum circu- 
lar velocity at the time of accretion, Vmax.acc 30kms~^, is 
abandoned. The outcomes presented in this section should 
not be considered as results per se. They are shown to point 
out the impact of the 30 kms~^ restriction. In this case, 
all satellites which can be detected at present, i.e. all sub- 
structures which have more than 100 particles at z=0, con- 
tribute to the averaging process. If a sub-halo falls below the 
100 particle limit we assume it as tidally dissolved. Since 
many satellites have low masses at the time of accretion 
tidal dissolution is more prominent in the unrestricted sam- 
ple. This is the key difference between the satellite sample 
with and without restriction. In the following we discuss the 
individual panels in Figure |4] and compare them to the cor- 
responding panels in Figure (2] 

a) The average accretion redshift, Zacc, is reduced for the un- 
restricted sample. This is expected since small mass satel- 
lites have shorter survival times. Thus, a large fraction of 
early accreted low mass systems is missing. However, the 
dependence on environment persists; 

b) The average number of satellites increases by a factor of 
< 2. This is a conseque nce of the shape of the (sub-)halo 
mass function (Bovlan- Kolchin et al. 2009); 

c) The average fraction between current maximum circular 
velocity and that at the time of accretion, Vmax/Vmax,acc is 
increased. This behaviour reflects the low average accretion 
redshift for the unrestricted satellite sample; 



6 Faltenbacher 



0.9 



O.E 



0.7 



0.6 



0.5 



a) 



10 



11 



,12 



lO^'- 10 
M,...[h-'MJ 



high density 
low density 
13 



10 



14 





3.0 




2.5 




2.0 


s 






1.5 






O 






1.0 




0.5 




0.0 




10 



b) 














hi0h Hpnsityr 




low density 


11 


jgl2 jq13 



14 



Mj...[h-'MJ 



0.95 
0.90 : 
^ 0.85 - 

> 

0.80 ; 
0.75 _ 



c) 



10 



1 1 



,12 



lO^'- 10 
M„[h-'MJ 



high density 
low density 
13 



10 



14 



0.2 
0.1 
0.0 
-0.1 
-0.2 
10 



d 



11 



high density 
low density 



10^^ 10 
M^..i[h-'MJ 



13 



10 



14 



2.0 



1.5 - \ 



1.0 - - 



high density 
low density 




0.5 [ 
10 



11 



1 ? 

10 10 
M„[h-'MJ 



13 



10 



14 



0.2 
0.1 
«i 0.0 
-0.1 
-0.2 
10 




1 1 



high density 
low density 



10^^ 10 
M„[h-'MJ 



13 



10 



14 



Figure 4. Same quantities as in Figure[2]but based on all satellites, i.e. without any restriction on the maximum circular velocity at the 
time of accretion. 



d) The average radial velocities, Wrad, are negative for 
all host masses. This is a consequence of the lower mass 
li mit of 100 particles for the sat ellites. As demonstrated 
in iFaltenbacher fc Mathewd l|2007l ) small mass satellites are 
most likely to fall below the resolution limit at their peri- 
centre passage causing an excess of negative radial velocities 
or inward moving satellites; 

e) The average velocity dispersions, ctsd, cTrad and atan, in- 
crease by a few per cent, independent of mass and environ- 
ment. In addition there is a substantial rise in the velocity 
dispersions of satellites in low mass hosts at high densities. 
The former effect can be explained by the lower average ac- 
cretion redshifts. More specifically, the increase of velocity 
dispersion is a result of the tidal dissolution of earlier ac- 
creted satellites which have on average lower velocities. After 
their dissolution the slow moving satellites do not contribute 
to the overall dispersion anymore which causes a positive 
bias of the velocity dispersion. This effect is more noticeable 
for the satellite sample without restriction of Vmax.acc since 
there tidal dissolution is more prominent. The increase of 
the velocity dispersion in low mass systems at high density 
is presumably an enhanced contribution of low mass inter- 
lopers in high density environments; 

f) The velocity anisotropy parameter, /?, decreases by ~ 
0.05, i.e. the tangential velocity component of the satellites 
is slightly increased for the sample without restriction on 
Knax.acc- lu addition, we find a distinct upward turn for 
low mass hosts in dense environments. The former is due 
to the fact that satellites on more radial orbits penetrate 
deeper into the potential well of the host. Consequently, they 
are more strongly exposed to tidal forces and get dissolved 
faster. The remaining satellite population is biased towards a 
more prominent tangential velocity component. The upward 
trend for low mass hosts at high densities may be explained 



by interlopers. Dynamically, they are not strongly coupled 
with the host. Assuming random motions of the interlopers 
one expects /? = which, in deed, is nearly approached. 

The comparison above highlights the fact that satel- 
lites are modelled by dark matter particles of a finite mass 
(6.9 X 10® /i~^Mq) and consequently can not be traced 
further or are lost when they get tidally striped below 
a given 'trustworthy' mass. Here we set this mass to be 
6.9 X 10* /i~^Mq corresponding to 100 particles. Satellites 
which get tidally striped below this limit 'disappear', i.e., 
they do not contribute to the average values (like accretion 
redshift, velocity dispersion etc.) for the satellite populations 
anymore. This causes a bias towards later accreted satellites 
and the consecutive effects. In particular the large number 
of approaching satellites with masses only slightly above the 
mass cut contributes to that bias. However, if the restriction 
Knax,acc ^ 30kms~^ is Imposed these effects are not appar- 
ent. In that case the satellites comprise on average a few 
1000 particles at the time of accretion. Even substantial tidal 
stripping does not push the satellite below the lower particle 
limit, i.e. all arriving satellites have extended survival times. 
Therefore, the results presented in Figure [2] describe phys- 
ical phenomena and are not caused by the finite masses of 
the simulation particles. 



4 CONCLUSION 

Taking advantage of the superior resolution of the 
Millennium-II Simulation we determine the impact of 
assembly bias on galaxy systems. We conclude with a 
recapitulation of the main results: 

1) The average accretion redshift depends on the host mass 
and environment. For host masses of 10^^ /i'^Mq we find 
Zacc ~ 0.7 which increases for host halos above 10^"^ h~^yiQy 
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to Zacc ~ 0.8. The average accretion redshift of satellites 
halos in high density environments is larger, Az ~ 0.1, 
compared to that of their counterparts in low density 
regions. The dependence of satellite accretion times on 
environment is most likely caused by the more tangentially 
biased velocity structure (see point 3) of satellites in high 
density regions which makes them more resistant to the 
hosts' tidal field. In contrast satellites in low density envi- 
ronments show more radially biased velocity structures. As 
a consequence total tidal dissolution of this satellites may 
happen faster causing a lower average accretion redshift of 
the surviving satellite population. 

2) Host halos above 10^^ /i^^M© show the same velocity 
dispersion as their satellite populations. This remains valid 
for lower mass halos in low density environments. However, 
the dispersion of satellites within low mass hosts at high 
densities exceed that of the hosts by ~ 30%. Which may be 
interpreted as contamination by high velocity interlopers. 

3) The velocity anisotropy of satellite halos is correlated 
with the mass of the host halo. More massive hosts show 
more radially biased satellite velocities. In addition we 
find a strong dependence of the velocity anisotropy on 
environment. Satellite velocities of host halos residing in 
high density environments tend do be less radially biased 
than those in low density environments. 

Our findings show that the internal dynamical structure 
of satellite systems is correlated with environment. We be- 
lieve that our approach allows to make similar predictions for 
real satellite galaxy systems. Thus the dynamical structure 
of satellite galaxy systems should depend on environment. 
Most notably satellites of host galaxies in less dense environ- 
ments show a more radially biased velocity structure. This 
behaviour may have implications for the accretion process 
of the host galaxies. 
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